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Abstract. The recent theoretical developments in the field of neutrino-nucleus interactions in the few-GeV region are 
reviewed based on the presentations made at the Nulnt09 Workshop. The topics of electron scattering and its connections with 
neutrino interactions, neutrino induced quasielastic scattering and pion production (coherent and incoherent) are covered, with 
special emphasis on the challenges that arise in the comparison with new experimental data. 
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INTRODUCTION 

Neutrino interactions with nuclei have received a consid- 
erable attention in recent years stimulated by the needs 
of neutrino oscillation experiments. A variety of the- 
oretical calculations have been performed for the dif- 
ferent reaction channels. At the same time, new high 
quality data are becoming available from MiniBooNE, 
MINOS, NOMAD and SciBooNE, and more is ex- 
pected from MINERvA, an experiment fully dedicated 
to cross section (CS) measurements. These problems 
have been addressed in the Sixth International Workshop 
on Neutrino-Nucleus Interactions in the Few-GeV Re- 
gion (Nulnt09) [1]. In this article, I present some of the 
theoretical highlights of this meeting, focusing on the 
most relevant reaction channels for accelerator-based ex- 
periments: quasielastic (QE) scattering and pion produc- 
tion (7rP); electron-nucleus (eA) interactions and their 
relevance for neutrino-nucleus (vA) scattering are also 
discussed. To begin, the motivations for doing theoreti- 
cal and experimental research on vA CS are reviewed. 

MOTIVATION 

Oscillation experiments J2l- Nowadays, the main rea- 
son for CS studies is the demand from oscillation ex- 
periments. Next generation ones aim at a precise de- 
termination of mass-squared difference Am^ and mix- 
ing angle 623 in Vj, disappearance measurements. The 
ability to reconstruct the neutrino energy is crucial for 
this program. Indeed, the oscillation probability P(v^ — > 
V T ) = sin 2 2023 sin 2 [Am\ 2 L/ (2E V )] depends on the neu- 
trino energy E v , which is not known for broad fluxes. The 
neutrino energy is usually reconstructed from charged- 
current quasielastic scattering (CCQE) events n — > 
pL ~ p, dominant at E v ~ 1 Ge V. If the energy and angle of 
the final lepton are measured, E v can be determined with 



two-body kinematics, but this procedure is exact only for 
free neutrons. As detectors are composed of nuclei, the 
reconstructed energy is smeared due to the momentum 
distribution of the bound nucleons. Moreover, the E v de- 
termination could be wrong for a fraction of events that 
are not CCQE ones but look identical to them in the de- 
tector. These are mainly A(1232) excitation events where 
the A is absorbed before decaying. Rejecting them re- 
quires a good command of final state interactions (FSI). 

High sensitivity searches for — > V e appearance, as- 
sociated with 0i3 and CP violation, also rely on good 
CS knowledge. The main background is neutral current 
(NC) 71° production because the electromagnetic show- 
ers from pion decay photons can be misidentified as elec- 
trons from v e induced CCQE v e n — > <?~ p if one of the 
photons is not identified. 

Astrophysics [3]. Neutrinos play an important role 
in astrophysical phenomena and carry information about 
the emitting sources. More precise measurements of the 
CNO neutrinos will provide a test for energy generation 
in stars and a better understanding of solar metalicity. 
The dynamics of core-collapse supernovae is controlled 
by neutrino interactions. The neutron rich environment 
of supernovae is a candidate site for r-process nucleosyn- 
thesis because radiated neutrinos convert neutrons into 
protons. These questions require a good knowledge of 
low energy neutrino production and detection CS. Al- 
though some neutrino CS of astrophysical interest could 
be investigated with beta beams, most of them cannot be 
measured and one relies on shell model or random phase 
approximation (RPA) calculations. 

Physics beyond standard model JHH]. Neutrino CS 
might be used to set bounds on nonstandard neutrino 
interactions. For example, deviations from universality 
in the Zw vertex could be accessed in deep inelastic 
scattering experiments at TeV energies. At low energies, 
a measurement of weak magnetism in v e p —> e + n may 
provide another test of vector current conservation. 



Hadronic physics. With high intensity neutrino 
beams it is possible to investigate the axial structure 
of the nucleon and baryon resonances. This informa- 
tion shall enlarge our view of hadron structure beyond 
what is presently known about electromagnetic form 
factors (FF) from JLab. MINERvA will probe the 
four-momentum transfer squared (Q 2 ) dependence of 
the nucleon axial FF with unprecedented accuracy and 
also study neutrino induced 7rP, which is dominated by 
resonance excitation. Another fundamental question is 
the strangeness content of the nucleon spin which can be 
best unraveled with NCQE v p{n) — > V p(n) reactions. 

Nuclear physics. Modern neutrino experiments are 
performed with nuclear targets. For nuclear physics this 
represents a challenge and an opportunity. A challenge 
because the above-mentioned prospects for more precise 
knowledge of neutrino and baryon properties require 
that nuclear effects are under control. An opportunity 
because vA CS incorporate a richer information than eA 
ones providing an excellent testing ground for nuclear 
structure, many-body mechanisms and reaction models. 



ELECTRON SCATTERING 

There are many theoretical methods that can be applied 
to vA interactions depending on the kinematic region: 
from shell model at low energies, through Fermi gas 
with hadronic degrees of freedom at intermediate ener- 
gies, to perturbative QCD at the highest ones. All of them 
have been applied to eA scattering. Moreover, some have 
been developed with the aim of understanding the large 
amount of good quality data acquired in decades of ex- 
perimental studies. A good description of electron scat- 
tering data is mandatory for neutrino interaction models. 

In addition, electron scattering data on proton and 
deuteron are used to extract electromagnetic nucleon 
elastic and transition (N — A, N — N*) FF, an input for the 
weak (V — A) hadronic currents since the V FF can be 
expressed in terms of the electromagnetic ones assuming 
isospin symmetry. Thanks to JLab data it has been estab- 
lished that proton electric and magnetic FF do not have 
the same behavior at large Q 2 [5]. These developments 
are incorporated in the FF parametrizations provided by 
Bodek et al. |g]. For resonance excitation, a unitary iso- 
bar model (MAID) has been used to extract the transition 
helicity amplitudes from the world data on pion photo- 
and electro-production for all four star resonances with 
masses below 2 GeV [7]. The helicity amplitudes can be 
mapped into the electromagnetic FF. The analysis reveals 
that the N — A(1232) transition is not purely magnetic. 
This has some impact on neutrino induced 7tP CS H . 

Let us consider now inclusive eA scattering, where 
only the final electron is detected. The simplest model 
for this reaction in the QE region, and the one imple- 



mented by default in all event generators used in the anal- 
ysis of neutrino data, is the relativistic global Fermi gas 
(FG) model. It assumes that the interaction takes place 
on single nucleons whose contributions are summed in- 
coherently [impulse approximation (IA)]. The struck nu- 
cleons have momentum distributions characterized by a 
Fermi momentum pf, and a constant binding energy Eb- 
Outgoing nucleons cannot go into occupied states (Pauli 
blocking) [9]. With two parameters (pf,£b), the main 
features of the inclusive QE CS can be explained. How- 
ever, a more careful view reveals the shortcomes of this 
simple picture. In particular, the global FG model over- 
estimates the longitudinal response. Nuclear dynamics 
must be taken into account flOll . 

The presence of nucleon-nucleon (NN) interactions 
implies that nucleon propagators are dressed with com- 
plex selfenergies E. In other words, the nucleons do not 
have a well defined dispersion relation but become broad 
states characterized by spectral functions 
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for both struck (holes) and outgoing (particles) nucleons. 
Sh includes an 80-90% contibution from single-particle 
states while the rest of the nucleons participate in NN 
interactions (correlations) and are located at high mo- 
mentum (see Fig. 4 of Ref. ITlll ). S p includes the ef- 
fect of the interaction of the outgoing nucleon with the 
medium, often accounted for with an optical potential 
(OP) within the Glauber approximation 111 ill . An alter- 
native approach [8] adopts a local FG model (pf (r) = 
[(3/2)^ 2 p(r)] 1 / 3 with p(r) the nuclear density), which 
introduces space-momentum correlations absent in the 
global FG (see Fig. 6 of Ref. (Htl). Instead of the constant 
binding, all nucleons are exposed to a density and mo- 
mentum dependent mean field potential. S p is considered 
in full glory, employing the low density approximation to 
calculate ImE, but it is argued that the correlated part of 
Si, plays a minor role in the description of inclusive CS 
and is neglected. With both frameworks a good descrip- 
tion of inclusive data in the QE region is achieved (see 
Figs. 2,3 of Ref. JH and Figs. 9,10 of Ref. H) improv- 
ing the result of the global FG. Furthermore, they can 
be extended to include resonances and nonresonant nP. 
While the agreement with data at the second (A) peak 
is also good, the description of the dip region between 
the QE and A peaks requires 2-particle-2-hole (2p — 2h) 
excitations from meson exchange currents (MEC) lfl3ll . 
The dip region is important for vA experiments because 
many CCQE-like events originate there. 

There is another class of IA relativistic models origi- 
nally developed for QE electron scattering and later ex- 
tended to neutrino scattering ffl [l5ll . The initial nu- 
cleons are treated as single-particle bound states whose 



wave functions are solutions of the Dirac equation with 
a (7-0) mean field potential. The various treatments of 
the final state include plane-wave IA, where the in- 
teraction of the outgoing nucleon with the medium is 
neglected, the so called relativistic mean field model 
(RMF), where the scattering wave functions are calcu- 
lated with the same energy independent real potential 
used for the bound nucleons, and distorted-wave IA with 
complex OP (DWIA); in this case the scattering states 
are obtained by solving the Dirac equation or using the 
Glauber model. DWIA models are successful in describ- 
ing a large amount of exclusive proton knockout (e,e'p) 
data but are not appropriate for inclusive scattering. The 
imaginary part of the OP produces an absorption and a 
reduction of the CS (see Fig. 1 of Ref. [ 14]), which ac- 
counts for the flux lost towards other channels. This is 
correct for an exclusive reaction but not for an inclusive 
one where all channels contribute and the total flux must 
be conserved. This approximation that retains only the 
real part for inclusive processes conserves the flux but is 
conceptually wrong because the OP has to be complex 
owing to the presence of inelastic channels ||l4Tl. An al- 
ternative is the Green function approach UoA where the 
imaginary part of the OP is responsible for the flux redis- 
tribution among different channels. This method guaran- 
tees a consistent treatment of both exclusive and inclu- 
sive reactions. An excellent description of the QE peak 
at different energies is achieved (Fig. 8 of Ref. IU6I1 ) al- 
though the transverse response is overestimated due of 
the lack of more complicated mechanisms such as MEC. 

Inclusive electron scattering data exhibit interesting 
systematics that can be used to predict vA CS. When 
the experimental (e,e') differential CS are divided by the 
corresponding single nucleon CS and multiplied by the 
global Fermi momentum, the resulting function 
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is found to depend on energy and 3-momentum trans- 
fers (co, \q\) through a particular combination, the scal- 
ing variable y/, and to be largely independent of the spe- 
cific nucleus (superscaling) IU7I1 . Scaling violations re- 
side mainly in the transverse channel. Therefore, an ex- 
perimental scaling function /( V 7 ) could be extracted by 
fitting the data for the longitudinal response. The experi- 
mental /( y/) has an asymmetric shape with a tail at pos- 
itive y/ (large co). The requirement of a realistic descrip- 
tion of the scaling function is a constraint for nuclear 
models. The relativistic FG model with exact superscal- 
ing gives a wrong symmetric shape for /(y/) while the 
RMF model reproduces it well. With the superscaling ap- 
proximation (SuS A) a good representation of the nuclear 
response can be obtained by embedding nuclear effects 
in the scaling function: the observables can be calculated 



with the simple relativistic FG model followed by the re- 
placement fpc — * fexp- The same strategy can be used 
to predict vA CS, minimizing the model dependence of 
the results. SuSA predicts 15 % smaller total CCQE CS 
compared to the relativistic FG. (Fig. 3 of Ref. IU8IP . It 
should be remembered that scaling fails at CO < 40 MeV 
and \q\ < 400 MeV due to collective effects. 



NEUTRINO INDUCED QE SCATTERING 

The CCQE scattering amplitude on a single nucleon is 
proportional to the product of the leptonic and hadronic 
currents. The hadronic one is given in terms of vector and 
axial FF Fa,p- Fp can be related to Fa using PCAC. Fa 

is usually parametrized as F A {Q 2 ) = gA (l + fi 2 /-^!) 2 \ 
with g A = 1.27 from j3 decay. For Ma, the world av- 
erage value from early experiments is Ma = 1 .026 ± 
0.021 GeV jl9tl and the one extracted from threshold 
% electroproduction data M A = 1 .069 ± 0.016 GeV d]. 
While the result of early neutrino experiments might be 
questionable due to the low statistics and poor knowl- 
edge of the neutrino flux, n electroproduction offers a 
solid indication that, at least at low Q 2 , Ma ~ 1 GeV. 

Modern experiments have started to provide a wealth 
of data on neutrino induced QE scattering for differ- 
ent energies and nuclear targets. MiniBooNE, running 
at (E v ) — 750 MeV on a CH 2 target, has collected 
the largest sample available so far for low energy V^, 
CCQE T20ll . After subtracting the non CCQE background 
using NUANCE [21], the CCQE data set was analyzed 
with the relativistic global FG model. The shape of the 
muon angular and energy distributions averaged over the 
Vu flux could be described with standard values of pf 
and Eb, but restricting the phase space for the final pro- 
ton by means of an ad hoc parameter K = 1 .019 ± 0.01 1, 
and taking M A = 1 .23 ± 0.20 GeV fH . This value of M A 
is considerably higher than the world average and the 
NOMAD result at high energies (3-100 GeV), also on 
12 C [M A = l.05±0m(stat)±0.06(syst) GeV] JH. A 
recent MiniBooNE reanalysis, using CC single nP data 
to adjust the simulation employed to subtract the back- 
ground, obtains k = 1.007 ±0.007 and an even higher 
M A = 1.35 ±0.17 GeV iB. 

The introduction of K offers a convenient way of 
parametrizing the low Q 2 reduction shown by the data 
within a simple FG model but its physical meaning is 
obscure. As explained above, the high Ma is also hard to 
understand. An alternative is that the observed increase 
of the CS at higher Q 2 might reflect the underlying nu- 
clear physics rather than the nucleon's. A more realistic 
description using spectral functions explains neither the 
low Q 2 reduction nor the high Q 2 increase l24l l25ll . In 
Ref. 12511 . it has been shown that the shape of the Q 2 



distribution extracted from data by MiniBooNE could 
be reproduced fairly well with Ma = 1 GeV by taking 
into account the renormalization of the electroweak cou- 
plings caused by the presence of strongly interacting nu- 
cleons (long range RPA correlations). This phenomenon, 
known as quenching, is well established in nuclear j3 de- 
cay and has proved crucial for a simultaneous descrip- 
tion of muon capture on 12 C and the low energy LSND 
CCQE measurement |26]. The problem is that quench- 
ing causes a reduction of the integrated CS already with 
respect to the FG result with = 1 GeV. The model 
prediction (a) — 3.2 x 10~ 38 cm 2 is much smaller than 
MiniBooNE's result of 5.65 x 10~ 38 cm 2 with an error 
of 10.8 % [23]. Other approaches like SUSA and RMF 
also find a reduction in the integrated CS. A promising 
solution to this puzzle has been proposed in a recent 
article i27l) . where the additional strength measured by 
MiniBooNE is explained by the contribution of 2/? — 2/z 
states that are not experimentally distinguishable from 
the standard \p—\h CCQE. More accurate comparisons 
to data (preferably inclusive) are still required: although 
the authors of Ref. 12111 stress that the dominant contri- 
butions to the 2p — 2h CS do not reduce to a modification 
of the A width, the background subtracted from the Mini- 
BooNE data includes a supression of % production that, 
in principle, accounts for pionless resonance decay ll2~lll . 

PION PRODUCTION 

Pion production in nuclei can be incoherent if the final 
nucleus is excited, or coherent if the nucleus remains in 
its ground state. The first step towards a good descrip- 
tion of %P on nuclear targets is a realistic model of the 
elementary reaction (on nucleons). In the few-GeV re- 
gion, the A(1232) excitation, followed by A — > nN is 
the dominant mechanism. Little is known about the ax- 
ial N - A FF, often denoted as Cj_ 6 . At moderate Q 2 , C5 
is the relevant one. The value of C£(0) can be related to 
the ANn coupling: c£(0) = gmnh/W^m N ) ~ 1 -2. The 
Q 2 behavior can be obtained by comparison to deuterium 
data from ANL and BNL bubble chamber experiments, 
with good statistics but large systematic errors due to the 
poor knowledge of the neutrino flux. In addition, there 
are nonresonant contributions that, close to threshold, are 
fully determined by chiral symmetry 12811 . The inclusion 
of these nonresonant terms led the authors of Ref. l28Tl 
to reduce C^(0) = 0.867 ± 0.075 to describe the data, 
but the fit was done to ANL data alone, which are sys- 
tematically below BNL ones. A combined analysis of 
both ANL and BNL data taking into account normaliza- 
tion uncertainties and deuteron effects, but not the non- 
resonant background, obtains C£(0) = 1.19±0.08 iB . 
Another theoretical description is based on a dynami- 



cal model of photo-, electro- and weak nP lf30h . Starting 
from an effective Hamiltonian with A — A couplings ob- 
tained with the constituent quark model (~ 30 % below 
the measured ones), the T matrix is obtained by solving 
the Lippmann-Schwinger equation in coupled channels. 
In this way the bare couplings get renormalized by me- 
son clouds. The predicted CS are in good agreement with 
data (Figs. 5-8 of Ref. JH). 

When 7TP takes place inside the nucleus, the elemen- 
tary cross section is modified. The most important in- 
medium change in this case is the modification of the 
A, whose mass gets shifted and its width increased due 
to absorption processes, mainly AN — > AW. The pro- 
duced pion interacts strongly with the nuclear environ- 
ment: it can be absorbed or scatter with the nucleons with 
and without charge exchange. At intermediate energies, 
a large number of states can be excited, so the description 
of the final system, in particular pion and nucleon propa- 
gation, requires a semiclassical treatment. The Giessen 
Boltzmann-Uehling-Uhlenberg (GiBUU) model allows 
to study these FSI in a realistic manner using transport 
theory in coupled channels. The effect of FSI on pion 
spectra appears to be large, especially for heavy targets 
(see for instance Figs. 13 and 14 of Ref. [31]). The ratio 
o(CC\n + ) / o(CCQE) as a function of the neutrino en- 
ergy has also been studied and compared to the observed 
(i.e. without FSI corrections) MiniBooNE experimental 
result [32]. The calculation should include all CCQE- 
like events in the denominator. It is found that the theory 
clearly underestimates the data at E v > 1 GeV (see Fig. 3 
of Ref. 13311 ). A similar result is obtained in Ref. [34] 
with a model that accounts for in-medium modifications 
of the elementary CS and propagates the pions through 
the medium with a cascade. Reference 112711 finds a good 
agreement for the observed ratio (bottom panel of Fig. 
19) but without including pion FSI. The latter will reduce 
the numerator and increase the denominator spoiling the 
agreement to some extent. 

Coherent nP occurs at very low Q 2 where the nucleus 
is less likely to break. It has a very small CS compared to 
the incoherent process, but relatively larger than coher- 
ent nP induced by photons or electrons due to the non 
vanishing contribution of the axial current at the relevant 
kinematics [35]. This reaction has attracted the attention 
of theoreticians because the low energy experiments find 
CS smaller than predicted by the pioneering model of 
Rein and Sehgal (RS) Q- The RS model uses PCAC in 
the Q 2 = limit to relate neutrino induced coherent kP to 
pion -nucleus (7tA) elastic scattering, which is modeled in 
terms of the pion-nucleon CS. By taking the Q 2 = limit, 
the RS model neglects important angular dependence at 
low energies 13711 . This, together with the fact that the de- 
scription of the tzA elastic CS is not realistic [37], results 
in CS well above the experimental data. As the energy in- 
creases, these deficiencies become less relevant and the 



RS model is revealed as a simple and elegant method. 
An alternative approach based on PCAC directly uses 
the experimental 7tA elastic CS IL381 l39ll . In this way the 
treatment of the outgoing pion is improved, but a spuri- 
ous initial pion distortion, present in %A elastic scattering 
but not in coherent nP, is introduced. With this method 
smaller and more compatible with the experiments CS 
are obtained. Microscopic approaches meant to work in 
the A region (and only there) have also been developed. 
They combine the A excitation picture of weak nP on the 
nucleon or the more complete models of Refs. IT281 l30ll 
with the A-hole model. As the nucleus remains in its 
ground state, a quantum treatment of pion distortion is 
feasible by means of the eikonal approximation l4fj|l . the 
Klein-Gordon l4ll l35ll or the Lippmann-Schwinger l42ll 
equations with realistic OR It is found that medium ef- 
fects and pion distortion reduce considerably the CS and 
shift the peak to lower pion-momenta (see for example 
Fig. 2 of Ref. 113510 . Nonlocalities in the A propagation 
were found relevant 14311 and have been explicitly taken 
into account in Ref. [42] but not in Refs. EJ l35ll . 
where they might have been partially included via the 
empirical A mass shift 14211 . The coherent nP CS is very 
sensitive to the value of (0) ll4Tll . A preliminary com- 
parison to MiniBooNE NCft data [44] indicates that a 
Cf(0) ~ 1.2 is preferred. 

CONCLUSIONS 

This is an excellent time for theoretical studies on vA 
interactions as new high quality data have started to ap- 
pear. A good understanding of (semi)inclusive vA (to- 
gether with eA) CS is required for the (model dependent) 
separation of different mechanisms: only then more pre- 
cise determinations of E v and of the NCn° background 
will be possible. The theoretical progress should also find 
its way to the simulations employed by the experiments. 
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